# Load data
nsduh <- read.csv("../Data/Modified_2022_NSDUH.csv", stringsAsFactors = TRUE)

# Subset youth only
youth <- nsduh[nsduh$is_minor, ] # subset for all minors
rm(nsduh) # dump original df

# Run OLS Models
padcigs.rolsmod <- lm(risk.pad_cigs ~
                        youth.parents_drugeduc +
                        youth.outclass_drugeduc +
                        youth.inclass_drugeduc,
                      data = youth)
bingedrink_daily.rolsmod <- lm(risk.bingedrink_daily ~
                                 youth.parents_drugeduc +
                                 youth.outclass_drugeduc +
                                 youth.inclass_drugeduc,
                               data = youth)
heroin_try.rolsmod <- lm(risk.heroin_try ~
                           youth.parents_drugeduc +
                           youth.outclass_drugeduc +
                           youth.inclass_drugeduc,
                         data = youth)
cocaine_monthlyuse.rolsmod <- lm(risk.cocaine_monthlyuse ~
                                   youth.parents_drugeduc +
                                   youth.outclass_drugeduc +
                                   youth.inclass_drugeduc,
                                 data = youth)


# Export intervention ordinal logits models
save(padcigs.rolsmod,
     bingedrink_daily.rolsmod,
     heroin_try.rolsmod,
     cocaine_monthlyuse.rolsmod,
     file = "../Results/Analysis_3.rdata")